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I .  INTRODUCTION 

Recent  advances  in  armor  technology  have  made  it  desirable  to 
examine  the  effectiveness  of  long  rod,  kinetic  energy  penetrators. 

As  the  length  to  diameter  ratio  of  the  rod  increases,  the  possibility 
arises  that  the  strong  axial  forces,  imposed  on  the  rod  by  impact 
on  a  target,  may  give  rise  to  unstable,  transverse  motions.  If 
tne  axial  load  is  imposed  for  a  long  enough  period  of  time,  the 
transverse  displacements  may  become  sufficiently  large  so  that  the 
rod  may  be  said  to  have  buckled.  Clearly,  severe  transverse  bend¬ 
ing  may  be  expected  to  degrade  the  penetration  process,  particularly 
in  the  case  of  multiple  impacts. 

It  is  the  purpose  of  this  report  to  estimate  the  conditions 
required  to  produce  transverse  instability  in  long  rod  penetrators. 
Penetration  mechanics  has  always  been  a  difficult  subject  for 
analysis,  and  frequently,  progress  has  followed  most  rapidly  from 
an  experimental,  rather  than  an  analytical,  approach.  In  this 
case,  however,  a  simple,  linear  analysis,  published  in  1966 
by  G.  R.  Abrahamson  and  J.  N.  Goodier  *,  may  be  used  to  make  explicit 
estimates  of  lateral  stability.  It  should  be  remarked  at  the  out¬ 
set,  however,  that  the  range  of  validity  of  the  theory  has  not  been 
completely  delineated,  since  Abrahamson  and  Goodier  reported  only 
a  limited  set  of  experiments  and  obtained  only  qualitative  agree¬ 
ment  in  some  cases. 

In  the  next  section,  the  theory  of  Abrahamson  and  Goodier  is 
summarized.  In  Section  III,  it  is  shown  how  to  apply  the  results 
of  Section  II  to  particular  design  cases.  In  Section  IV,  the 
direction  of  future  work  is  indicated. 


II.  SUMMARY  OF  THEORY** 

When  a  rod,  moving  with  high  axial  velocity,  strikes  a  target, 
large  axial  stresses  may  be  induced  at  the  impact  end.  If  the  rod 
could  respond  elastically,  a  discontinuous  increase  in  stress, 
given  by  equation  (1),  would  prop  agate  down  the  bar  with  speed  c  . 

A  o  =-p  cg  A  u  (1) 


♦References  are  listed  on  Page  (21). 
**6ee  reference  1  for  complete  details. 
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Here,  A  a  is  the  increase  in  stress,  p  is  the  density  of  material 
in  the  rod,  cg  is  the  elastic  bar  speed  (ce  =  /E ^/p  ,  where  E^ 

is  the  elastic  modulus),  and  A  u  is  the  instantaneous  change  in  velocity 
of  the  end  of  the  rod.  To  see  how  large  A  a  may  be  in  a  representative 
case,  consider  a  steel  rod  striking  a  steel  plate  with  an  impact  velocity 
of  4000  fps.  In  an  elastic  impact  of  materials  with  identical  properties, 
one  might  estimate  A  u  to  be  2000  fps.*  With  E  taken  as 

30  x  10  psi,  and  the  specific  gravity  taken  as  7.8,  the  calculated 
increase  in  stress  comes  out  to  be  about  3.5  x  10^  psi,  which  is 
approximately  sixteen  times  the  yield  stress  of  high  strength  tool  stools 


(e.g.,  bearcat  steel  o 


2.2  x  105 


psi). 


Clearly,  the  impact  end  of  the  rod  cannot  sustain  such  stresses  for 
long,  and  certainly,  such  a  high  stress  level  cannot  propagate  into  the 
body  of  the  rod.  Instead,  the  end  of  the  rod  is  worn  away,  a  stress 


equal  to  a 


the  yield  stress,  will  propagate  down  the  bar  with 


a  speed  c  ,  the  elastic  bar  speed,  and  stress  levels  above  a  will 
e  y 

begin  to  propagate  down  the  bar  with  a  speed  dictated  by  the  slope  of 

the  strain  hardening  curve.  The  speed,  ,  of  these  additional 

stresses  is  often  quite  small  compared  to  the  elastic  bar  speed; 


c  /c  =  0.13  to  0. 19 
P/  e 


in  the  cases  considered  by  Abrahamson  and  Goodier.  Furthermore,  bending 
that  occurs  after  passage  of  the  elastic  wave  will  depend  on  the  strain 
hardening  modulus,  E^  .  For  the  purpose  of  estimating  stability,  with 

the  preceding  discussion  as  motivation,  the  following  simplified  model 
is  constructed.  Refer  to  Figure  1. 

Assume  that  the  stress-strain  curve  of  the  bar  is  bilinear  with 

elastic  modulus  E  ,  yield  stress  o  ,  and  strain  hardening  modulus 
e  y 

E^  (see  Figure  2,  for  example).  Since  axial  loading  from  longitudinal 

wave  propagation  occurs  much  more  rapidly  than  transverse  motion,  it 
is  assumed  further  that  the  yield  stress  is  applied  continuously  to  the 
end  and  in  the  body  of  the  rod,  and  that  additional  bending  occurs  as 
a  small  perturbation  so  that  all  longitudinal  fibers  experience  increasing 
load.  Rotatory  inertia  and  shear  deformation  are  neglected,  and  simple 
beam  theory  is  applied.  The  bending  moment  M  is  proportional  to  the 
second  derivative  of  lateral  displacements  y 

M  =  -Eh  A  k2  &L  (2) 


*This  estimate  is  based  merely  on  the  instantaneous  velocity  after  impact 
of  like  materials  and  has  nothing  to  do  with  residual  velocity. 
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where  A  is  the  cross  sectional  area  of  the  rod  and  k  is  the  radius 
of  gyration.  Shear  force  Q  induces  lateral  acceleration. 


ia, 

3x 


P  A 


(3) 


Balance  of  moments  requires  that 


Q  *  t  (y  *  yJ 


m 

3x 


(4) 


where  P  is  the  axial  load, 
tion. 


and  y  is  the  initial  deflec- 

*  O 


In  the  actual  boundary  valve  problems  to  be  examined,  the  finite 
rod  is  assumed  to  be  sufficiently  long  that,  in  the  analysis,  it  may  be 
replaced  by  a  semi-infinite  rod.  Abrahamson  and  Goodier  assumed  that  both 
deflection  and  moment  vanish  at  x  =  0;  y  (0,t)  =  0,  M  (0,t)  =  0. 

These  quantities  must  vanish  at  infinity  as  well.  Two  initial  conditions 

were  used:  either  yo  (x)  =  0,  ^  o;  or  yo  (x)  t  0,  =  o. 


Presumably,  these  two  conditions  were  assumed  to  give  representative 
results.  In  any  case,  the  main  results  turn  out  to  be  estimates  of 
times  for  buckling  to  occur  and  characteristic  wavelengths  in  the 
buckled  pattern.  The  detailed  distribution  of  buckles  is  not  predicted. 

The  various  parameters  are  expressed  in  the  following  nondimensional 

way. 


w 


K 


t 

k 


2 

s  = 


(5) 

(6) 

(7) 


Nondimensional  displacements  ire  w  and  w0,  and  nondimensional 
distance  and  time  are  £  and  t  respectively.  In  terms  of  the 
variables  defined  by  equations  (5)  and  (6),  the  equation  of  motion 
is 


IV 


+  w"  +  w  =  -w„ 


(8) 
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w 


where  w'  =  3w/3  5  and  w  =  3w/3x.  Since  it  is  assumed  that  both 
deflection  and  moment  vanish  at  5  =  0,  the  solution  may  be  expressed 
as  a  Fourier  sine  transform 


w  (5,  t}= 


t) 


sin  n  5  d  n  • 


(9) 


where  g  (n,  x)  must  satisfy  the  ordinary  differential  equation 


g  -  n  (1  - 


n2) 


g  =  o 


(10) 


Initial  transverse  velocity  is  supposed  to  arise  "...  for  instance, 
from  elastic  bending  waves  following  a  slightly  oblique  impact",  or 
possibly  from  previous  impacts.  If  we  assume  that  w0(5)  =  0  and  that 
initial  conditions  are  given  by 


w  (5,  0)  =  0 


w  (5,  0) 

then  g  (5,  n)  is  given  by 


8  (n)  sin  n  5  d  n  , 


(ID 


g  (5  ,  n)  =  v0  c  (n)  x 


p~*  sinh  p 


-1 

p  sin  p 


,  if  n  <  1 


,  if  n  >  1 


(12) 


where  p  =  xn 


n2)1'2  . 


Abrahamson  and  Goodier  now  argue  that  large  displacements  occur 
when  g  is  large.  At  each  fixed  valve  of  x  ,  this  occurs  when 

p  is  large  and  p  always  has  a  maximum  at  n  =  n  =  — 

'/2 
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The  amplification  factor  of  0  (n)  in  equation  (12), 


sinh 


P 


2  1/2 
n 


sinh  t n 1 1  - 


f 


plotted  for  0  <_  n  ^  1  with  a  somewhat  arbitrary  vertical  scale,  shows 

no  strong  peak  at  t  =  3,  but  by  x  =  6,  a  peak  is  forming,  and  by 

x  =  9  the  peak  is  well  developed  (see  Figure  3). 

This  description  does  not  change  much  if  w0  (5)  t  0  with 
initial  conditions  taken  to  be  w  (5,  0)  =  $  (5,0)  =  0  .  For  this  case, 
the  preferred  value  of  n  is  not  independent  of  t  ,  but  decreases 
with  increasing  t  .  A  peak  in  the  amplification  factor  again  forms 
by  t  =  6  and  is  well  developed  by  x  =  9,  but  now  the  peak  occurs  at 

a  value  of  n  some  12%  to  15%  larger  than  before.  There  is  no 

theoretical  reason  for  singling  out  the  values  t  =  3,  6,  or  9,  but 
these  values  do  seem  to  be  useful  for  correlation  of  the  experimental 
results  in  reference  1. 

The  net  result  is  that  buckles  with  nondimens ional  wave  number 
n  may  be  expected  to  have  formed  by  x  =  6  and  to  be  fully  developed 
by  x  =  9,  provided  that  loading  has  been  maintained.  From  equations 
(5)  aid  (6),  this  means  that  buckles  of  half  wavelength  X  will  form 
at  times  t  given  by 


sn 


t  = 


where  x  lies  roughly  in  the  range  6  <_  x  <_  9,  and  buckling  will 
become  progressively  more  severe  for  still  larger  times. 


nk 


<2Eh/V 


1/2 


(13) 


kx 


2 

sc  a 

p  y 


(13) 


III.  DESIGN  APPLICATIONS 

It  must  be  recognized  that  the  conclusions  of  the  last  section 
rest  on  an  incomplete  asymptotic  analysis.  The  magnitude  and  location 
of  buckles  cannot  be  estimated  without  an  explicit  statement  of 
initial  values  and  more  extensive  analysis.  For  design  purposes. 
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however,  the  results  of  the  last  section  may  be  sufficient  if  they  are 
interpreted  conservatively. 

For  instability  to  show  up,  the  loading  time  must  be  greater  than 
the  time  to  buckle.  There  are  two  unloading  mechanisms: 

1.  The  rod  penetrates  and  passes  through  the  target  plate,  or 
else  it  is  stopped.  In  either  case,  the  load  is  released  after  time 
tp  ,  which  may  be  estimated  by 


Here  h  is  the  target  thickness  for  complete  penetration  or  the 
depth  of  penetration  otherwise,  and  vm  is  the  mean  projectile  speed 

during  penetration.  If  uniform  deceleration  is  assumed,  vm  is  the 

average  of  the  striking  and  residual  velocities.  The  nondimen sional 
time  to  penetrate  is 


TP 


(14) 


Comparison  of  xp  with  any  other  time  t 


shows  that  the  ratio 


is  >,  =,  or  <  1 


according  as 


or 


>1  / 

1  Eh 


h 

F 


(15) 


2.  The  elastic  wave  reflects  from  the  rear  end  of  the  rod  where 
the  stress  drops  back  to  zero.  The  unloading  wave  then  propagates 
at  the  elastic  bar  speed  back  to  the  impact  end.  The  time  required 
for  the  round  trip  is 


t 


R 


2L 

c 

e 


2L 
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The  nondimensional  time  is 


Tn  = 


X.  it 


R  t  k 
h  e 


(16) 


Comparison  of  td  with  any  other  time  t  shows  that  the  ratio 
R 


according  as 


or  <  1 


or 


T_ 

2 


E 

e 


a 

y 


(17) 


If  loading  is  maintained  at  the  impact  end,  the  loading  wave  will  be 
reinitiated.  It  seems  likely  that  the  resultant  cyclic  loading  and 
unloading  for  a  given  length  of  time  will  be  less  likely  to  cause 
buckling  than  sustained  loading  for  the  same  length  of  time.  Further¬ 
more,  if  buckles  do  form,  they  will  probably  be  less  severe.  These 
facts  should  be  accounted  for  in  establishing  design  rules. 

Equations  (15)  and  (17)  are  the  primary  results  with  design  implica¬ 
tions  in  this  memorandum.  For  t  =  3,  6,  and  9,  these  two  equations 
are  plotted  in  Figures  4  and  5.  Radius  of  gryation  k  has  been 
replaced  by  the  symbol  d/4.  For  a  solid  rod  of  uniform  material 
throughout,  d  is  the  diameter. 


Figures  4  and  5  are  to  be  used  as  follows.  For  a  given  long  rod 

penetrator  and  a  given  target,  the  numbers  1,  d  (or  k) ,  E.  ,  E  , 

n  g 

0^  ,  and  p  are  known.  E^,  E£  ,  and  ay  should  be  obtained,  preferably, 

from  carefully  run  compression  tests.  Then  c^  =  ^Ef/P  may  be  computed, 

and  for  a  given  impact  velocity,  the  numbers  vm  and  h  must  be 

estimated.  Next,  the  ordinates  and  abscissas  for  the  two  figures 
may  be  computed  and  a  point  located  in  each  figure.  The  range 
of  values  of  Xp  or  x^  that  corresponds  to  each  wedge  shaped 


region  has  been  indicated  on  the  figures.  Examination  of  the  experi¬ 
mental  results  of  Abrahamson  and  Goodier  show  that  tentative  inter¬ 
pretation  of  the  points,  so  located  in  the  various  wedges  of 
Figures  4  and  S  ,  might  be  as  follows.  In  each  case,  Xp  is  considered 
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first,  aiul  interpretation  is  then  conditioned  by  Tp  .  The  design 
rule  is  stated  first  in  italics,  followed  by  a  few  explanatory  comments. 

1)  If  x  p  <  2,  then  the  rod  in  probably  safe  for  any  xp  .In 

other  words,  the  time  of  penetration  is  sufficiently  short  that  the 
impact  end  unloads  before  buckles  have  time  to  develop  no  matter  how 
long  the  rod. 

ii)  If  2  <  Tp  <  6",  then  the  rod  is  probably  safe  if  xp  <  9, 

but  it  nay  be  only  marginally  safe  if  x  >  9.  The  load  is  now 

maintained  at  the  impact  end  for  a  long  enough  time  that  buckling  may  be 
incipient  unless  unloading  by  reflection  from  the  free  end  also  occurs. 

If  buckles  do  begin  to  form,  the  amplitudes  (maximum  deflection  from  the 
straight  rod)  will  be  small.  In  subsequent  impacts,  the  incipient 
buckles  will  tend  to  be  amplified. 

iii)  If  0  <  ip  <  9 ,  then  the  rod  is  probably  safe  if  t  <  6 , 

and  marginally  safe  for  d  <  xp  <  9.  If  xp  >  9 ,  the  rod  may  still  be 

marginally  safe.  The  same  comments  as  for  ii)  apply  here  as  well,  but 
now  since  time  of  penetration  is  more  severe,  the  time  for  reflection 
judged  to  be  safe  has  been  reduced. 

iv)  If  xp  >  9 ,  then  the  rod  may  be  safe  for  x p  <  «J,  but  for 

2  <  Xp  <  6',  it  is  only  marginally  safe.  For  6  <  xp  <  9  the  rod  may 

still  be  marginally  safe  provided  xp  <  2  xpi  but  if  xp  >  2  xp 

or  if  xp  >  9 ,  then  buckles  will  probably  form.  The  longer  the  loading 

time  Tp  or  the  larger  the  ratio  x 

buckling  will  be.  In  this  case,  the  load  is  maintained  at  the  impact 
end  for  a  sufficiently  long  time  to  develop  buckling.  For  short 
rods,  however,  the  load  is  periodically  released  from  the  rear  end  so 
that  the  transverse  motion  must  be  reinitiated  each  time  the  axial 
loading  passes  by.  As  tp  and  Tp  increase,  there  is  enough  time 

for  buckles  to  develop  either  during  one  passage  of  the  axial  wave 
(Tp  >9,  xR  >  9)  or  during  multiple  passages  of  the  axial  wave 

(Tp  >  9,  6  <  Xp  <  9,  but  Tp/Tp  >  2  or  3) .  There  is  nothing  in  the 

analysis  that  indicates  the  conditional  effect  of  Xp  for  large 

Tp  (i.e.  Tp  >  9)  as  stated  here.  Rather,  it  is  inferred  from  examination 

of  the  experimental  results  of  Abrahamson  and  Goodier,  and  it  is  con¬ 
sistent  with  those  results. 

The  interpretation  of  the  last  four  paragraphs  is  only  tentative 
because  the  analysis  of  Abrahamson  and  Goodier  does  not  predict  the 
severity  of  buckling.  Furthermore,  their  experiments  were  not  designed 


p^Ap  ,  the  more  severe  the 
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to  find  criteria  for  the  marginal  onset  of  buckling,  but  rather  were 
intended  to  demonstrate  the  wavelength  relationship  once  buckling  had 
definitely  occurred,  finally,  their  rods  were  made  from  6061  -  T6  and 
2024  -  T3  aluminum,  rather  than  high  strength  steels  or  high  density 
materials,  and  it  has  not  been  demonstrated  that  the  theory  will  then 
apply.  Therefore,  the  above  conclusions  cannot  be  drawn  directly  from 
the  data,  but  represent  an  extrapolation  of  some  sort. 

Still  further  extrapolation  is  required  to  apply  these  results 
to  technologically  important  cases.  Figure  2  shows  an  idealized  stress- 
strain  curve  for  bearcat  steel.  The  ratio 


c 

P/ 


'c 

e 


0.152  , 


computed  from  the  curve,  falls  in  the  same  range  as  in  the  experiments 
of  Abrahamson  and  Goodier  (0.13  to  0.19)  so  that  the  rate  of  plastic 
loading  compared  to  elastic  loading  is  similar.  On  the  other  hand, 
the  important  ratio  a  =  0.319  for  bearcat,  but 

only  0.13  to  0.25  for  aluminum  in  reference  1,  so  that  the  stress 
strain  curves  are  only  approximately  similar.  Furthermore,  impact 
conditions  may  vary  widely  from  those  used  in  reference  1.  As  noted 
in  the  introduction,  the  ratio  A  a/a  -  16  for  an  impact  velocity  of 


4000  fps  in  bearcat,  but  the  maximum  in  all  the  experiments  of 
Abrahamson  and  Goodier  is  A  a/a  -  8  ,  so  that  the  excess  energy 


available  above  the  elastic  strain  energy  absorbed  at  yield  is  greater 
for  typical  impacts  with  bearcat  rods.  In  addition,  the  ratio 
v  /Cp  for  bearcat  (cp  -  2560  fps)  may  well  be  in  the  range  of  1.5  or 

2.0,  whereas  the  maximum  in  reference  1  is  only  0.29.  Finally,  in 
the  experiments  there  were  no  significant  penetrations,  and  the  impact 
ends  of  the  rods,  although  slightly  mushroomed  in  some  cases,  were 
never  actually  worn  away. 


In  spite  of  these  reservations,  figures  4  and  5,  if  used  conser¬ 
vatively,  should  be  helpful  at  least  for  preliminary  design. In  figure  4, 
the  vertical  and  horizontal  lines  have  been  drawn  for  bearcat  steel 
with  ratios  of  h/d  =  1,  2,  5,  10,  and  15  and  mean  penetration  velocities 
of  v  =  2000,  4000,  and  6000  ft. /sec.  In  figure  5,  vertical  lines 
corresponding  to  bearcat  steel,  6061-T6  aluminum,  and  2024-T4  aluminum 
are  shown.  Other  materials  could  be  plotted  as  well  of  course. 


IV.  FUTURt  WORK 

At  the  end  of  the  last  section,  it  was  indicated  that,  in  many 
respects,  applications  may  well  fall  outside  the  limits  of  experimental 
verification.  Further  experimentation  can  readily  show  whether  or  not 
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the  scaling  is  shown  to  be  valid,  however,  there  remains  the  fact  that 
it  applies  only  to  the  case  of  normal  incidence  without  yaw.  These 
two  factors,  oblique  incidence  and  yaw,  will  certainly  introduce  trans¬ 
verse  forces  on  the  impact  end  of  the  rod  and  may  introduce  a  transient 
moment  as  well.  Axial  compression  will  tend  to  amplify  the  bending 
in  these  cases.  An  investigation  of  the  effect  of  transverse  end 
forces  or  end  moments  on  solutions  of  equation  (8)  is  currently  underwa)  and 
will  be  reported  upon  in  the  near  future. 
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Figure  5.  Time  of  Load  Release  from  Reflections  at  Rear  of  Projectile 
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